Goos-Hanchen and Imbert-Fedorov shifts of a nondiffracting Bessel 
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Goos-Hanchen and Imbert-Fedorov shifts are diffractive corrections to geometrical optics that have been 
extensively studied for a Gaussian beam that is reflected or transmitted by a dielectric interface. Propagating 
in free space before and after reflection or transmission, such a Gaussian beam spreads due to diffraction. We 
address here the question how the Goos-Hanchen and Imbert-Fedorov shifts behave for a "nondiffracting" 
Bessel beam. © 2010 Optical Society of America 
OCIS codes: 240.3695, 260.5430. 



It has been known since a long time that the behavior 
of a finite-diameter light beam in reflection and transmis- 
sion at a dielectric interface differs from the predictions 
of geometrical optics. Due to diffractive corrections the 
beam is shifted in directions parallel and perpendicular 
to the plane of incidence pQ . The parallel shift is known 
as the Goos-Hanchen (GH) effect [2] and the transverse 
shift as the Imbert-Fedorov (IF) effect [3J4] . These effects 
have been extensively studied [5HZ| , not only for total in- 
ternal reflection which is the context wherein the GH and 
IF effects were originally addressed, but also in partial 
dielectric reflection and transmission [5]. We note that 
the IF effect is closely related to the Spin Hall Effect of 
Light (SHEL) [5], [5UTUUTT] . Further generalizations con- 
cern angular varieties of the GH and IF effects; these are 
observed in the far field of the reflected (or transmitted) 
beam [TJj . Recently, also the influence of Orbital Angu- 
lar Momentum of the incident beam on these diffractive 
shifts has been investigated [T3"rfT5] : for the transverse 
case this may be called the Orbital Hall Effect of Light 
(OHEL) [16] . 

The diffractive origin of these effects raises the ques- 
tion how they behave when the incident beam is a 
so-called nondiffracting Bessel beam. Such beams were 
perceived by Durnin et al. as propagation-invariant so- 
lutions of the free-space scalar wave equation (TT1IT5] . 
These solutions have amplitudes proportional to Bessel 
functions. The zero-order Bessel beam has a bright cen- 
tral maximum ( "needle beam" ) which propagates in free 
space without diffr actional spreading; the higher-order 
beams have a dark central core. Most work on Bessel 
beams is restricted to the paraxial limit [19l[20] but also 
the nonparaxial case (described by the Helmholtz wave 
equation) has been reported [2TH22] . 



Ideal Bessel beams have infinite transverse diameter 
and can therefore not be generated experimentally. How- 
ever, there exist several experimental methods to gen- 
erate finite-diameter approximations to a Bessel beam; 
these propagate over a finite axial distance in a non- 
diffracting manner (i.e. over distances much larger than 
the Raleigh length corresponding to the needle-beam di- 
ameter) [23] ■ So, one may speculate that such a needle 
beam corresponds to a geometrical-optics ray that would 
not show GH and IF shifts. Verification or rebuttal of 
this speculation requires proper theory; this is reported 
in the present paper. 

Let us begin by briefly recalling to the reader what 
a Bessel beam is. The scalar m th -order Bessel beam is 
a cylindrically symmetric monochromatic optical beam 
whose electric field has the following form: 



E(R,<p t z) = J m {K R)e imip e lz ^ k °- K % 

= A(R,(py z V k2 o- K Z, 



(i) 



where m is an integer number that fixes the value of 
the orbital angular momentum (OAM) of the beam, and 
(i?, tp, z) are the cylindrical spatial coordinates defined 
with respect to the main axis of propagation z: 



x = Rcosip, 
y = R sin tp. 



(2) 



For a needle beam one has m = 0. In Eq. ([T]) ko > 
and < Kq < ko are two independent parameters, 
where Kq determines the angular width §0 of the central 
lobe (cone) of the corresponding Bessel function via the 
definition 



K = fc sin#o, 



(0 < ??o < tt/2). 



(3) 
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Throughout this Letter we will consider only paraxial 
Bessel beams characterized by the condition 



sin$ = K /k < 1. 



(4) 



The physical meaning of the angle #o is illustrated in Fig. 
1 below. It should be noticed that while E(R, tp, z) is an 




Fig. 1. Illustrating the physical meaning of the angle in- 
exact solution to the Helmholtz equation (d 2 + d 2 + d 2 + 
k^)E = in free space, the amplitude A(R,<p) satisfies 
the reduced equation (d 2 + d 2 + Kq)A = 0. 

For actual calculations of both GH and IF shifts it is 
convenient to work in Fourier space and calculate the 
Fourier transform A(K, cf>) of the amplitude A(R, tp) = 
J m {K R)e lm * as: 

J m (K R)e imifl = i- J A(k x ,k v )e iK R dk x dky, (5) 



where 



A{k x ,k y ) = A(K, 



^S(K-K )e" n t (6) 
r'«o 



with K = (k x + ky) L/z , K R = xk x + yk y = 
KRcos((f> — tp), and k x = K cos cf>, k y = Ksmcj). 
It is worth noticing that in the literature Eq. ([5]) is 
often written in spherical coordinates (ko, •&,(/)) with 
K = fcosintf, K S(K-K ) = S (i? - i? ) /costf , and 
dk x dk y = k 2 sin i? cos i? d$d0. 

Having written explicitly the Fourier representation of 
a scalar Bessel beam, we can now proceed as in |16j and 
write the Fourier amplitude of a vector Bessel beam as: 

A(k X ,ky) -> A(k X , ky) = f ± (K)A(k X ,ky), (7) 

where f±(K) = f — k(k- f), with fe = k x x + k y y + (k 2 — 
K 2 y/ 2 z and / = f p x+f s y, \f\ 2 = 1. Here, according to 
[To] , the three unit vectors {x, y, z} form a right-handed 
Cartesian reference frame attached to the incident beam 
propagating along the axis z. 

At this point the GH and IF shifts for a Bessel beam 
impinging at the angle 8 upon a planar interface may be 
straightforwardly calculated by using the formulas given 



in Sec. Ill of Ref. [24]. Once again, following Ref. [16] we 
define the "intrinsic" (namely beam-independent) lon- 
gitudinal and transverse shifts as, respectively, 



. d In r a 

' 88 



.R\ 



(8) 



and 



Y '-f.V 



where r\ = R\ exp(i<j)\), A <G {p, s}, with the prime 
indicating derivatives with respect to the incidence angle 
8. Moreover, we define the relative reflected energies w p 
and if, as: 



„ 2 p 2 



a p 2 Rp 2 + a s 2 R s 2 



, w s 



. 2p 2 

a s U s 



a p 2 R p 2 + a s 2 R s 2 



, (10) 



and the complex-valued longitudinal and transversal 
shifts S and $ respectively, as: 



S = WpX p + w s X s , 



* = w p Y p + w 3 Y 3 , (11) 



where 

„ ,r^s _ ap 2 R v 2 4>p + a-s 2 R s 2 4 l s' 
tie (z.) — 2 2 , 

a p z R p + a s z R s 

-j- — ,^ dp RpRp ~\~ ®>s R-s-Rs 



s 2 R~7~ 

2 , r, 2\ 



a p 2 R p + a 2 R 



Re(tf) 



a p a s cot 8 (R p + R s ) sin 77 

a p 2 R p 2 + a 2 R s 2 
a p a s cot 8 [2R p R s sm(rj — 4> p + 4> s )] 



a 2 R 



a s 2 R s 



, , a a» cot 8 (R D 2 — R s 2 ) cos 77 
Im (*) = - p 



a s 2 R s 



(12a) 
(12b) 

, (12c) 
(12d) 



From Ref. [2SJ it immediately follows that the "tra- 
ditional" GH and IF shifts evaluated for a fundamen- 
tal Gaussian beam of waist wq and Raylcigh range 
L = /coiwq/2 can be expressed in terms of the formulas 
given above as: 

k Q (x r ) = Re(S) + ^lm(S), (13a) 
k (y r )=Re(y) + jlm(¥), ( 13b ) 

where the real parts of 5 and VP give the spatial shifts of 
the beam, and the imaginary parts furnish the angular 
GH and IF shifts defined as d(x r )/dz r and d(y r )/dz r , 
respectively. It should be noticed that here and in the 
subsequent formulas, according to Ref. [16], the three 
Cartesian coordinates {x r ,y r , z r } are referred to a ref- 
erence frame attached to the reflected beam of central 
wavevector fco, with z r directed along fco- X straightfor- 
ward calculation shows that in the case of a m th -order 
Bessel beam, Eqs. ([T3]) become: 



k o(x r )\ Zr=0 = Re(H) - mlm(*), (14a) 
*b(»r>L =0 = Re(*)+mIm(S), (14b) 
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for the spatial part, and 

^=sin^Im(S), (15a) 

oz r 

= sintfglm(tf), (15b) 

oz r 

for the angular part. These formulas are the main re- 
sult of this Letter. Before proceeding with the discussion 
of these formulas, a caveat is in order here. They were 
derived in straightforward manner by analogy with the 
Gaussian beam case. However, while a Gaussian beam is 
describable by means of normalizable functions, a Bcsscl 
beam does not. Luckily, in the practical calculation of 
GH and IF shifts, infinities present in the first-order mo- 
ments of the electric field energy density distribution are 
exactly (and, unambiguously) compensated by the in- 
finities given by the electric field energy density inte- 
grate over the whole space. Thus, the non-normalizable 
nature of (theoretical) Bessel beams does not represent 
a problem. 

Some relevant issues follow from Eqs. (|14II15[) above. 
First, since for a fundamental Gaussian beam of an- 
gular aperture 6q one has l/(koL) = Oq/2, then from 
Eqs. (fTS"]) with sin^o ~ i?o, it follows that the angular 
shift of a Bessel beam is about twice the correspond- 
ing shift of a Gaussian beam. Second, Eq. (fbfj) shows a 
spatial/angular mixing analogous to the one present for 
a Laguerre-Gauss beam of OAM m [15] . Third, in par- 
ticular for the case m = 0, the formulas derived above 
have profound consequences for the (nondiffracting) core 
of the Bessel beam, i.e. the needle beam. The key point 
is that in first order perturbation theory (with respect 
to the expansion parameter $o) the (full) Bessel beam is 
not deformed upon reflection, so it translates rigidly as a 
whole. Since, as we saw above, the full Bessel beam shows 
the standard GH and IF shifts this must also be the 
case for its central nondiffracting core. The only escape 
from this conclusion is to leave the paraxial approxima- 
tion and to go to second-order (and higher) perturbation 
orders (i.e. relatively strong focusing); this reduces the 
length over which diffraction is effectively absent to a 
propagation length of the order of the diameter of the 
full Bessel beam. However, a treatment of beam shifts in 
this regime is outside the scope of this paper. 

A. A. acknowledges support from the Alexander von 
Humboldt Foundation. 
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